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Abstract. We study standing waves for a model of nonlinear Schrodinger equation on a graph Q. 
The graph is obtained joining N halflines at a vertex, i.e. it is a star graph. At the vertex an interaction 
occurs described by a boundary condition of delta type with strength a ^ (which includes the free 
or Kirchhoff case a = 0). The nonlinearity is of focusing power type. The dynamics is given by an 
equation of the form i4r^! t = H^ t — |^ t | 2 ' t1 $ r t , where H is the Hamiltonian operator which generates 
the linear Schrodinger dynamics on the graph. We show the existence of several families of standing 
waves, i.e. solutions of the form $ t = e^*^ where ^ is a vector amplitude parametrized by a 
frequency u>; they exist for every sign of the coupling at the vertex, attractive [a < 0) and repulsive 
(a > 0) and for every u> > ■ The number of such families depends on the number N of the edges 
of the graph. Furthermore, we determine the ground states, which are the standing waves whose 
amplitude minimizes the action on the natural (Nehari) constraint, and order the various families of 
standing waves according to their increasing action, thus determining a nonlinear spectrum of the 
problem. Finally, we show that the ground states are orbitally stable for every allowed ui if the 

2 

nonlinearity is subcritical or critical; for supercritical nonlinearities there exists an u* > p such that 

2 

the ground states are orbitally stable for < us < lo* . The minimization of the action requires 
technical subtleties. In particular, we develop a non-trivial adaptation of decreasing rearrangements 
on a star graph, and obtain a generalization of the Polya-Szego inequality for the kinetic energy of a 
rearranged state. 
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1. Introduction 

In the present paper a rigorous analysis of the stationary behaviour of nonlinear Schrodinger equation 
(NLS) on a graph is given, beginning from the simplest type of unbounded graph, the star graph. In 
a previous paper j2] the authors studied the behaviour in time of an asymptotically solitary solution 
of NLS resident on a single edge of the graph in the far past, and impinging on the vertex with various 
types of couplings, giving a quantitative analysis of reflection and transmission of the solitary wave 
after the collision at the junction. Here we concentrate on a different phenomenon, namely the 
existence of persistent nonlinear bound states on the graph (localized, or pinned nonlinear modes), 
and on their orbital stability, when an attractive interaction is present at the vertex. 
Let us briefly give a collocation of the model in the physical context. Generally speaking, one can 
consider the NLS as a paradigm for the behaviour of nonlinear dispersive equations, but it is also an 
ubiquitous model appearing in several concrete physical situations. The main fields of application 
which we have in mind are the propagation of electromagnetic pulses in nonlinear media (typically 
laser beams in Kerr media or signal propagation in optical fibers), and dynamics of Bose-Einstein 
condensates (BEC). We are interested in the way solutions of NLS are affected by the presence of 
inhomogeneities of various type. The propagation on the line in the presence of defects has been a 
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subject of intense study in the last years and it gives rise to quite interesting phenomena, such as 
defect induced modes [13j [TTJ [7], i.e. standing solutions strongly localized around the defect. The 
presence of defect modes affects propagation by allowing trapping of wave packets, as experimentally 
shown in the case of local photonic potentials in [28J. On the other hand, nonlinearity can induce 
escaping of solitons from confining potentials, as demonstrated in [3T]. A last interesting phenomenon 
is the strong alteration of tunneling through potential barriers in the presence of nonlinear defocusing 
optical media |34j . In this paper we consider NLS propagation through junctions in networks. 
For example, when the dynamics of a BEC takes place in essentially one-dimensional substrates 
("cigar shaped" condensates) or a laser pulse propagates in optical fibers and thin waveguides, the 
question arises of the effect of a ramified junction on propagation and on the possible generation of 
stable bound states. The analysis of the behavior of NLS on networks is not yet a fully developed 
subject, but it is currently growing. Concerning situations of direct physical interest we mention 
the analysis of scattering at Y junctions ("beam splitters") and other network configurations ("ring 
interferometers") for one dimensional Bose liquids discussed in (33]. Some more results are known 
for the discrete chain NLS model (DNLS), see in connection with the present paper the analysis in 
|2"9] . Other recent developments are in [TSJ E2]- In particular, in the paper [TS] scattering from a 
complex network sustaining nonlinear Schrodinger dynamics is studied in relation to characterization 
of quantum chaos. 

With these phenomenological and analytical premises in mind we would like to construct a mathe- 
matical model capable to represent, in a schematic but rigorous way, the propagation and stationary 
behaviour of a nonlinear Schrodinger field at a junction of a network. We begin to give the needed 
preliminaries to rigorously define our model. We recall that the linear Schrodinger equation on graphs 
has been for a long time a very developed subject due to its applications in quantum chemistry, nan- 
otechnologies and more generally mesoscopic physics. Standard references are [TQl [TU |2H 1251 I2"5] . 
where more extensive treatments are given. Here we recall only the definitions needed to have 
a self-contained exposition. We consider a graph Q constituted by N infinite half lines attached 
to a common vertex. The natural Hilbert space where to pose a Schrodinger dynamics is then 
L 2 (Q) = 0^ 1 L 2 (IR + ). Elements in L 2 {Q) will be represented as function vectors with components 
in L 2 (IR + ), namely 




We denote the elements of L 2 (Q) by capital Greek letters, while functions in L 2 (K + ) are denoted 
by lowercase Greek letters. We say that \& is symmetric if ip^ does not depends on k. The norm of 
L 2 -functions on Q is naturally defined by 

N 

i=i 

From now on for the L 2 -norm on the graph we drop the subscript and simply write || • ||. Accordingly, 
we denote by (•, •) the scalar product in L 2 {Q). 
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= max 

00 l^j^N 



y j||L°°(R+)- 



Analogously, given 1 ^ r ^ 00, we define the space L r (Q) as the set of functions on the graph whose 
components are elements of the space f/(R + ), and the norm is correspondingly defined by 

N 

i=l 

Besides, we need to introduce the spaces 

N N 

H\Q) = 0# X (R + ) H\Q) = 0ff 2 (R+), 

j=i 3=1 

equipped with the norms 

N N 

\Wm = \Wh* = 52Wi\\h<vi+y (1-1) 

i=l i=l 

Whenever a functional norm refers to a function defined on the graph, we omit the symbol Q. 
When an element of L 2 {Q) evolves in time, we use in notation the subscript t: for instance, ^> t - 
Sometimes we shall write ^{t) in order to highlight the dependence on time, or whenever such a 
notation is more understandable. 

The dynamics we want to set on the graph is generated by a linear part and a nonlinear one. We 
begin by describing the linear part. 

Fixed a G R, we consider a Hamiltonian operator, denoted by H and called S graph or 5 vertex, 
defined on the domain T>(H) 

N 

V(H):={*eH 2 (g)s.t. Vi(0) = ... = Vjv(0), ^)^(O) = a^i(O)}, (1.2) 

8=1 

where ip[ denotes the derivative of the function ipi with respect to the space variable related to the 
i-th edge. The action of the operator H is given by 

h^> = ; 

V -k 

The Hamiltonian if is a selfadjoint operator on L 2 (Q) ([23J) and generalizes to the graph the ordinary 
Schrodinger operator with 5 potential of strength a on the line [9]. Similarly to that case, the 



interaction is encoded in the boundary condition. The case a = in (1.2) plays a distinguished role 
and it defines what is usually given the name of free or Kirchhoff boundary condition; we will indicate 
the corresponding operator as H°. Notice that for a graph with two edges, i.e. the line, continuity of 
wavefunction and its derivative for an element of D(H°) makes the interaction disappear; this fact 
justifies the name of free Hamiltonian. A 5 vertex with a < can be interpreted as the presence of 
a deep attractive potential well or attractive defect. This interpretation can be enforced by showing 
that, as in the case of the line, the operator if is a norm resolvent limit for e vanishing of a scaled 
Hamiltonian H e = H° + aV e , where V e = W{^) and V is a positive normalized potential on the 
graph (see pj] and reference therein). The attractive character shows in the fact that for every 

2 

a < a (single) bound state exists for the linear dynamics, with energy —^2. On the contrary, on 
a Kirchhoff vertex no bound states exist, the spectrum is purely absolutely continuous, but a zero 
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energy resonance appears. Finally we recall that in the case of repulsive delta interaction a > 0, 
which is however of minor interest here, there are not bound states nor zero energy resonances. 
The quadratic form E hn associated to H is defined on the finite energy space 

£ = V{E lm ) = {V e H\Q) s.t. Vi(0) = ■ ■ ■ = Vjv(O)} 

and is given by 

N "+<*> n . 1 _ n 



J] / |^(x)| 2 rfx+-|^(0)| 2 = -||^ 
1=1 ^° 
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The corresponding bilinear form is denoted by B(-, •) and explicitly given by 

fl(*,$) := -^(^',0/)l 2(K+ ) + |^(O)0i(O). (1.3) 

i=l 

As a particular case, the quadratic form E°' hn associated to H° is defined on the same space, that is 
V(E°' lin ) = £, and reads 

N ,+oo 



i . . r + °° i 

E 0,Un m= J2 / M(X)\ 2 dx = -{WW 2 
1 i=l ^0 ^ 



Now let us introduce the nonlinearity. To this end we define G = (Gi, ■ ■ ■ , G^) '■ C n — > C n where G 
acts "componentwise" as Gj(C) = fi'dODO f° r a suitable g : 1R + —> E and C = (C?) ^ C n . 
We are interested in the special but important case of a power nonlinearity of focusing type, so we 
choose g(z) = — |;?| 2Ai , n > . 

After this preparation it is well defined the NLS equation on the graph, 

= H%-\%\ 2 »% (1.4) 

where /i > 0. This abstract nonlinear Schrodinger equation amounts to a system of scalar NLS 
equations on the halfline, coupled through the boundary condition at the origin included in the 
domain (1.2). 



In Section 2 we show that for /i > well-posedness of the dynamics described by equation (1.4) 
(in weak form) for initial data in the finite energy space £ holds true. Moreover, if < fi < 2 then 
the solution exists for all times and blow-up does not occur. Finally, as in the standard NLS on the 
line, mass M(\&) = |||^|| 2 and energy E[$f] are conserved, where 

m = l\\n 2 -^\\n&i+%\Mo)\* 

and analogously, in the case a = 0, for the Kirchhoff energy 

-■OhT/l _ 1 lhTr'l|2 1 lhTrll 2 M+2 



E"m = -w\\ 2 —11*11^X2 • 

L J 2 II II 2^ + 2" ll2M+2 



After setting the model and its well posedness (see Section |2[), w e turn to the main subject of 
this paper, existence and properties of standing wave solutions to (1.4). Standing waves are solutions 
of the form 

%(x) = e wt * w (x) . 
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The function \[? u is the amplitude or the profile (with some abuse of interpretation) of the standing 
wave, and we will frequently refer to the set of ^ as to the stationary states of the problem. 
The amplitude ^ satisfies the stationary equation 

IH> W - |^| 2 ^ = -a;*, , tu>0. 

This equation has a variational structure. 
Let us define the action functional 

s u [v] = Em + = ^i*'n 2 + |n*r - ^^wnZXl + f i^iwi 2 . 

The Euler-Lagrange equation of the action is the stationary equation above. The action S u , defined 
on the form domain S of the operator H, is unbounded from below. Nevertheless, it is bounded 
on the so called natural (or Nehari) constraint G £ s.t. = 0}, where I u [^f] = \\^f'\\ 2 — 

ll^llKa + W II^H 2 + «l^i(0)| 2 - Note that J w (* w ) = S^(* w )* w , and thus the Nehari manifold is a 
codimension one constraint which contains all the solutions to the stationary equation. One of our 
main results is the following theorem. 

Theorem 1 (Existence of minimizers for the Action functional). Let fi > 0. There exists a* < 
such that for —Nyfu < a < a* the action functional S w constrained to the Nehari manifold admits 
an absolute minimum, i.e. a ^ ^ such that = and S u [^f] = inf {S^^] : = 0}. 

So the action admits a constrained minimum on the natural constraint for every uj > if the 
strength a of the 5 interaction at the vertex is negative and sufficiently strong. The strategy of the 
proof, which is a consequence of results of Section 3 and Section 4, makes use of non trivial elements 
and we give here some remarks. To get the existence of the minimum one has at a certain point to 
compare the action S 1 ^ with a < with the Kirchhoff action S®. 

In Section 3 we prove that the Kirchhoff action, while bounded from below on its natural 
constraint, has no minimum (see [3] for an analogous phenomenon affecting the constrained energy 
functional). As a matter of fact, the infimum can be exactly computed and it is achieved as the limit 
over a sequence of functions which escape at infinity on a single edge. A main step in establishing 
the previous picture and in applying it to the 5 case, is the exact calculation of the infimum and 
the identification of the minimizers of the free action; to this end one exploits an extension and 
generalization of the classical properties of symmetric rearrangements of LP and H 1 functions to the 
case of graphs. 

In Section 4 we prove Theorem 1. The analysis follows in part proofs of similar results for 
singular interactions on the line given in [T?l [6], with major modifications due to the fact that in 
this case the comparison with the free case is not standard. In particular the upper bound in a 
given in the statement of Theorem 1 is a consequence of the fact that one needs the condition 
infS^ < infS 1 ^ to guarantee the existence of an absolute minimum in the constrained action, and 
this to penalize situations analogous to escaping minimizers of the free action. A sufficiently strong 
attractive interaction at the vertex allows to satisfy the previous condition. 

We conjecture that the action has a local constrained minimum that is larger than the infimum when 
the condition on a fails, but presently we do not have a proof of this fact. 

In Section 5 an explicit construction of all the stationary states of the problem is obtained, by 

2 

solving the stationary equation for every value of a. It turns out that for every N and u > there 
exist families of stationary states of different action and energy, which can be ordered in j to 
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form a nonlinear spectrum (the family is unique only in the case N = 2, i.e. the line). The state of 
minimal action is the ground state ^ W} q, which is of course the solution to the constrained minimum 
problem for the action just discussed. The others are excited states, and they exists, for every j, 

2 

when uj > (jv° 2 j) 2 • ^ ee Theorem 4 for a complete description. 

Finally, in Section 6 we study the stability of ground states. Stability is an important requisite 
of a standing wave, because at a physical level unstable states are rapidly dominated by dispersion, 
drift or blow-up and so are undetectable (instability of NLS with a 5 potential on the line is studied, 
partly numerically, in [26]). The concept of stability, due to gauge or U(l) invariance of the action, 
is orbital stability. The solutions remain close to the orbit e* e \I/ Wi o of the ground state for all times 
if they start close enough to it. The framework in which we study orbital stability of the ground 
state is the mainstream of Weinstein and Grillakis-Shatah-Strauss theory, which applies to infinite 
dimensional Hamiltonian systems such as abstract NLS equation when a regular branch of standing 
waves u ty w<Q (not necessarily ground states) exists, which is our case. 

According to this theory, to guarantee orbital stability one needs to verify a set of spectral conditions 
on the linearization of the NLS around the ground state, and a slope (or Vakhitov-Kolokolov in the 
physical literature) condition concerning the behavior in u of the L 2 -norm of the ground state. Some 
adaptation of standard methods is needed to treat the singular character of the interaction at the 
vertex, but in fact it turns out that in the range over which it is valid the previous theorem on the 
variational characterization of the ground state ^^o, the spectral conditions and the slope condition 
are encountered for every u and every nonlinearity in the range /i G (0, 2]. 

Theorem 2 (Orbital stability of the ground state). Let \i G [0,2], a < a* < 0, u > jfa. Then the 
ground state ^f Wi o is orbitally stable in £ . 

The proof of this result is contained in Section 6. Notice that one has orbital stability of the 
ground state in a range of nonlinearities which includes the critical case. An analogous phenomenon 
occurs in the case of the line, previously treated in [26]. This marks a difference with the case of 
a free (a = 0) NLS on the line, where one has orbital instability in the critical case. Finally, the 
proof of the previous theorem (see Remark 6.1) shows that for supercritical nonlinearities \x > 2 

2 

the ground state is orbitally stable for not too large u: there exists a threshold u* > ^ such that 

2 

one has orbital stability for the ground state W Wi0 with < uj < uj* and orbital instability in the 
opposite case. 

Appendix A contains a theory of symmetric rearrangements on star graphs, which is a technical 
novelty. More precisely, the classical inequalities stating conservation of LP norms and domination 
of kinetic energy are proved. This last property, i.e. the Polya-Szego inequality, is particularly 
interesting because it changes with respect to the case of the line through the presence of a factor 
which takes into account the number of edges of the graph, and this fact is crucial in the previously 
described analysis of action minimization on a star graph. We stress the fact that the theory of 
rearrangements so constructed is a general tool and it is in principle applicable to more general or 
different problems. 

We end this introduction with a few open problems and future directions of study. Concerning 
technical issues, a different strategy from the one here pursued in the analysis of ground states 
and their stability is minimization of energy at constant mass (see the classical paper [T5] and for 
models related to the present one [8] ) ; it requires a non trivial extension to graphs of concentration- 
compactness method and will be studied in a forthcoming paper [T]. Nothing is known up to now 
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about stability properties of the branches of excited states which exist for every N > 2 and 
sufficiently high cu; this is a subject of special interest because there are only few cases where excited 
states of NLS equations are explicitly known. The authors plan to study this issue in a subsequent 
paper. Finally it would be interesting, and perhaps a difficult task, the extension of the analysis here 
given to different classes of graphs, possibly with non trivial topology. This will need the development 
of new technical tools, both concerning variational analysis and stability properties. 

Acknowledgments The authors would like to thank Reika Fukuizumi for several discussions. R. 
A. was partially supported by the PRIN2009 grant Critical Point Theory and Perturbative Methods 
for Nonlinear Differential Equations. 



2. Well-posedness of the model 

For our purposes it is sufficient to prove that the solution of the Schrodinger equation is uniquely 
defined in time in the energy domain and that energy and mass are conserved quantities. This section 



is devoted to the proof of these conservation laws and of the well-posedness of equation (1.4). In fact 



along the proofs we shall always work with the weak form of (1.4), namely 

q? t = e - iHt ty + i [ e~ iHit - s) \ty s \ 2 ^ s ds. (2.1) 
Jo 

We consider the problem of the well-posedness in the sense of, e.g., [H], i.e., we prove existence and 
uniqueness of the solution to equation (2.1 ) in the energy domain of the system. Such a domain turns 
out to coincide with the form domain of the linear part of equation ( |1.4 ). We follow the traditional 



line of proving first local well posedness, and then extending it to all times by means of a priori 
estimates provided by the conservation laws. Proceeding as in [2] where the cubic NLS is treated, 
we show the well-posedness of the dynamics for any /i > 0, i.e. local existence and uniqueness for 
initial data in the energy space. Moreover, we will prove that if < fi < 2, then the well-posedness 
is global, i.e. the solution exists for all times and no collapse occurs. For a more extended treatment 
of the analogous problem for a two-edge vertex (namely, the real line with a point interaction at the 
origin) see [5]. 



We endow the energy domain £ with the if 1 -norm defined in (1.1). Moreover we denote by £* the 
dual of £, i.e. the set of the continuous linear functionals on £. We denote the dual product of 
r G £ * and \l/ G £ by (r, In such a bracket we sometimes exchange the place of the factor in £* 
with the place of the factor in £: indeed, the duality product follows the same algebraic rules of the 
standard scalar product. 

As usual, one can extend the action of H to the space £, with values in £*, by 

(H^ 2 ) := B[#i,tf 2 ], 

where B[-,-] denotes the bilinear form associated to the selfadjoint operator H. 
Furthermore, for any \1/ G £ the identity 

^e~ iHt ^ = -iHe~ im ^ (2.2) 
dt y ' 
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holds in £* too. To prove it, one can first test the functional J^e lHt ty on an element H in the operator 
domain £, obtaining 



d 
db 



lim 



e iH{t+h)-z: _ e iHt-z: 

h 



(■V,iHe iHt Z) = {-iHe- lHt V,Z) 



-iHt,T, "\ 



Then, the result can be extended to S e £ by a density argument. 



Besides, by (2.2), the differential version (1.4) of the Schrodinger equation holds in £ 



In order to prove a well-posedness result we need to generalize standard one-dimensional Gagliardo- 
Nirenberg estimates to graphs, i.e. 



1^1 



C\\ty'\\*-p\\V\\s + ', (2.3) 

where the C > is a positive constant which depends on the index p only. The proof of ( 2.3[ ) 
follows immediately from the analogous estimates for functions of the real line, considering that any 
function in if 1 (IR + ) can be extended to an even function in if 1 (IR), and applying this reasoning to 
each component of \1/ (see also [30, 1.31]). 

Proposition 2.1 (Local well-posedness in £). 



Let > 0. For any \l/o G £, there exists T > such that the equation (2.1) has a unique solution 
^ e C°([0,T),£) nC l ([0,T),£*). Moreover, eq. fl2~T] ) has a maximal solution ^ m 
interval of the form [0, T*), and the following "blow-up alternative" holds: either T 



defined on an 
■ oo or 



lim 



\ w t \\£ 



+oo, 



where we denoted by \|>™ ax the function \]/ max evaluated at time t. 

Proof. We define the space X := L°°([0,T),£), endowed with the norm H^Ha" 
Given ^ e £, we define the map G : X — > X as 



su Pte[o,T) ||*tlU- 



G$ := e - iH ^ + i [ e- iH( ~ s) \$ s \ 2, *$ s ds. 
Jo 



We first notice that the nonlinearity preserves the space £. Then by 
using Holder and Gagliardo-Nirenberg inequalities, one obtains 



'I ^ C\(f)\^\(j)'\ and 



s\\S 



^ C||$ 



|2/j+1 



so that 



\\GQ\\ X ^\\*o\\£ + C [ \\^s\\T +1 ds ^ \\* \\e + CT\m% 

J o 

Analogously, given $,H G £, one has 

||G$ - GH|U < CT (||$||^ + ||H||^) ||$ - S|U 



2^+1 



(2.4) 
(2.5) 



We point out that the constant C appearing in (2.4) and (2.5) is independent of $, and 5. 
Now let us restrict the map G to elements $ such that || <^> || ^ ^ 2||\l/o||£- From (2.4) and (2.5), if T 



1 2m 



then G is a contraction of the ball in X of radius 



is chosen to be strictly less than (8C||^o llg 
2||^o|U; an d so ) by the contraction lemma, there exists a unique solution to (2.1) in the time interval 
[0, T). By a standard one-step boostrap argument one immediately has that the solution actually 
belongs to C°([0, T), £), and due to the validity of (1.4) in the space £* we immediately have that 



the solution * actually belongs to C°([0, T), £)) n C T [$,T),8 



The proof of the existence of a maximal solution is standard, while the blow-up alternative is a 
consequence of the fact that, whenever the £-norm of the solution is finite, it is possible to extend it 
for a further time by the same contraction argument. □ 

The next step consists in the proof of the conservation laws. 

Proposition 2.2 (Conservation laws). 

Let /i > 0. For any solution * G C°([0,T),S) n C 1 ([0,T),5*) to the problem ( |2l| >, the following 
conservation laws hold at any time t: 

M[%] =Af[* ], £[*t] = E[* ]. 

Proof. The conservation of the L 2 -norm can be immediately obtained by the validity of equation 
(1.4) in the space £: 



d 
dt 







by the self-adjointness of H. In order to prove the conservation of the energy, first we notice that 
(\Pt, H^f t ) is differentiable as a function of time. Indeed, 



[(% +h ,H% +h )-(%,H%)\ 



h 



h 



and then, passing to the limit h — > 
d 



-2Im(|* t | 2 "*„H* t ), 



where we used the self-adjointness of H and (1.4). Furthermore, 
d 



dt 



(* t ,|*t| 2 ^) 



d 
dt 



2(/i + l)Im(|^| 2 ^,#^>. 



(2.6) 



(2.7) 



From (2.6) and ( |2.7 ) one then obtains 



and the proposition is proved. 



□ 



Corollary 2.1 (Global well posedness). 

Let < /i < 2. For any G £> equation (2.1) /ias a unique solution ^> G C°([0,oo),£) fl 
^([0,00),^). 



Proof. By estimate (2.3) with p = oo and conservation of the L 2 -norm, there exists a constant C, 
that depends on \l/ only, such that 

E[V ] =E[V t ] > ^WtW* ~ CWT 

Therefore a uniform (in t) bound on \\^' t \\ is obtained. As a consequence, one has that no blow-up 
in finite time can occur, and therefore, by the blow-up alternative, the solution is global in time. □ 
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3. Variational Analysis: the Kirchhoff vertex 

In this section we compute the infimum of the action functional for the Kirchhoff case. As often 
in this framework, the action functional is unbounded from below and we have to restrict it to the 
Nehari manifold, or natural constraint manifold, in order to have a functional bounded from below. 
The knowledge of the infimum of the constrained action will be a key ingredient in the next section 
in the proof of the main theorem. 

The strategy of the computation of the infimum is standard: first we derive a lower bound and then 
we show that this lower bound is optimal by means of a minimizing sequence. In the derivation of 
the lower bound symmetric rearrangements are used. Using this technique we can map the initial 
variational problem into a variational problem with symmetric functions which can be reduced to a 
problem on the halfline providing the required estimate. 

The minimizing sequence shows in fact that the constrained action exhibits a sort of spontaneous 
symmetry breaking in the Kirchhoff case. That is, although the functional is symmetric, the mini- 
mizing sequence is localized on a single edge. 

As defined in the introduction, in the Kirchhoff case the action functional is given by 

s2M = E°m + coMm = i\\n 2 + f n*n 2 - ^^wnili (3.1) 

while the Nehari functional J°, reads 

m = \\n 2 -\\nzii+u\\n 2 - 

The Nehari manifold is defined by G £, \& 7^ s.t. = 0}. The action restricted to the 

Nehari manifold will be named reduced action and is given by 

m = asm - \iim = 2^ii*i$s ( 3 - 2 ) 

It is understood that the domain of all the functionals is always £. 
Theorem 3 (Infimum of the Action for the Kirchhoff case). 

The infimum of the action functional S® restricted to the Nehari manifold is given by: 

pi 

inf{SM s.t G £ , ^ 0, = 0} = d°(w) = (ji+ / (1 - t 2 )ht. (3.3) 

Jo 



Proof 

The proof of (3.3) is divided into two parts: first we derive a lower bound for S®, then we prove that 
the lower bound is optimal by means of a minimizing sequence. 

In order to derive a lower bound, we consider an auxiliary variational problem with symmetric 
functions. This is done by using the rearrangements on the graph which are discussed in Appendix 

\E 

Let $ G £ and let $* be its symmetric rearrangement. We known that <3>* is positive, symmetric and 
$* G £. Moreover, by Theorems [6] and A.l we have 

2 



1 2/^+2 



WW > 



(3.4) 
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Therefore for $ G 8 such that /°[$] = we have 

A||$*'|| 2 - ||$*||^ + W ||^||2 ^ jo^ = (3 _ 5) 

and 

£[$] = £[$*]. (3.6) 

Taking into account (3.2), and the above properties of $* one can enlarge the domain in the following 
way in order to lower the infimum, 



inf{S£[$] s.t. $ G S, $ ^ 0, J°[$] = 0} = inf{S[$] s.t. $ G £, $ ^ 0, J°[$] = 0} 

^ inf{£[$] s.t. $ G £, $ ^ 0, $ symmetric, ^||$'|| 2 - H^H^S + ^ll $ f < 0}. (3.7) 
Under the scaling, $(•) ~» A 1//2 $(A-), A > 0, the last variational problem scales as 

inf s.t. $ G £, $ ^ 0, $ symmetric , ^||$'|| 2 - H^HaJ+a + w ll $ ll 2 < °} = 

inf A M s.t. $ G £, $ ^ 0, $ symmetric , ^A 2 ^')! 2 - A^||$|| + w ll $ f < °| • 

It is convenient to choose A as 

4 /A/A ^ 

—A 2 = A- so that A = ^-J 

in order to reconstruct a Nehari manifold with a rescaled u; as constraint. Moreover due to the 
symmetry of $ we have 

(^J 2 " inf | s.t. $ G £, $ ^ 0, $ symmetric , ||$'|| 2 - H^lgja + w f^J 1 1 <^> 1 1 2 ^ 



2u 



JV f Vl 2 " inf | ^^||0|liC + 2 2(R+) s.t. G ^(R+) , 



7^ 0, H'\\ 2 L2{R+) - ||0||^ + 2 2(R+) + u l^pj " M < j . (3.8) 

It is convenient to introduce a variational problem on the half line and an auxiliary variational 
problem on the line. Let d (oj) and d lme (u) be defined in the following way: 

rf ha V) = mf {a^lMlSKw B -*- G ^ °> WW) " IMK+V) + fc#|&<* +) < o 

^ ine M = inf {^^||0|| 2 C + 2 2(K) s.t. G ^(K), 0/0, ||0'||| 2(R) - |H|g+ 2 a(Il) + w|M&(E) < 

Notice that the following inequality holds true: 

2rf half (w) ^ rf linc (w). (3.9) 



2 M 
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Indeed, by absurd, assume 2<i half (c<;) < d hne (u) and let <p n be a minimizing sequence for the problem 
on the halfline. We can extend <p n by parity and obtain a sequence <p n G H l (§L) such that 



n|li 2 ( 



ll^nllL 2 (R) — ll^n|ll,2M+2( R ) + ^ 

Passing to the limit one would obtain 

d lhic (u) > liminf 



sC 



1 2/^+2 
I L 2 ^+ 2 ( 



2\\<Pn" W 



L 2 M+ 2 (R+)- 



7 1 1 2/1+2 
J n\\ jJ.^+21 



Id 



half/ 



2/i + 2' 

which contradicts our absurd hypothesis. Therefore 1/2 d hnc (u) provides a lower bound for the 
variational problem we are interested in. On the other hand the exact expression of d lme (u) can be 
easily obtained from known results (see [H] Ch. VIII) , and it is given by: 

-l 



d linc (co) = (n + l)tu$ 



i i i i 

2 



[l-t 2 )»dt. 



(3.10) 



Taking into account (3.7), (3.8), (3.9) and (3.10) we can conclude 



2fi 

N f N\ 



2fj. 



V+2 
2/i 



'\-t 2 Y»dt 



1 M+ 2 
(/i+ljt'U 2 ^ 



'1 -t 2 )»dt. 



(3.11; 



Estimate (3.11) closes the first part of the proof. Now it is sufficient to exhibit a sequence of trial 

functions <£>„ satisfying the constraint and such that S[Q n ] — >■ (yU+1) ^(l— t 2 )»dt. We consider 

a sequence of soliton-like functions escaping to infinity, i.e. 

($n)t(aO = 



0„(x) = s (a; - n)x(x) i = 1 
i^l 



(3.12) 



where 4> s is defined in Appendix |B| by (B.l ) and x is a C 
for ^ x ^ 1 and x( x ) 



function such that 0^x^l?X = 



1 for x ^ 2. The sequence <3> n belongs to £ but does not satisfy the 



constraint /°[$ n ] = 0. It is straightforward to check that 



l$nl 



2^+2 



> C 



(3.13) 



where the r.h.s. of (3.13) depends on u and \i. In the remaining part of the proof we shall not make 
explicit the dependence on uj and fi of the constant appearing in estimates. Let 5 n be defined by 

Kl| 2 + ^ll$n|| 2 



i 



1$ ii 2m+2 



(3.14) 



It is straightforward to check that J°[5 n $ n ] = 0. Then in order to prove (3.3), it is sufficient to prove 
that 

-l 



Now we prove that 



lim S[5 n $ n ] = (/J, + l)iu;^. 

n— >oo 



lim <L 



:i - t 2 )»dt 



(3.15) 
(3.16) 
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We have 



iff. i|2/i+2\ 2 M 
1^*1112^+2 



1$ ii 2m+2 



and by (3.13), it is sufficient to prove that 

lim 11$: 'II 2 - 



w||$J 



1$ || 2 ^+ 2 - n 

I II 2^+2 _ u ' 



(3.17) 



(3.18) 



Taking into account (|3.12|) and (|B.2|) and integrating by parts one has 

ira a +wii*jr 



i$ ii 2 ^+ 2 i < 



ci \<p s (x -n)\ 2 \x"(x)\\x{x)\dx + c / \<p s {x - n)\4>' a {x - n)\\x (x)\\x(x)\dx = K n . (3.19) 

'0 JO 

The remainder TZ n can be estimated using the exponential decay of <p s and <fi' s in the following way 



x)\ 2 dx + c / \(j) s (x)(j)' s (x)\dx ^ c 

n—l Jn—1 



e cx dx ^ ce" 



(3.20) 



This proves (3.16) while ( 3.15[ ) is reduced to prove that 



Jim S[$ n ] = + / CI - /-')- ;; r//. 

'o 



The last equality follows by dominated convergence and (|B.5|): 



lim S[$ n ] = lim 



n-s>oo 2/i + 2 Jg 



5 (x — n)x(x)| 2/i+2 cix = lim 



2/i + 2 



|0 s (x)x(x + n)| 2/i+2 c/x 



2/i + 2 



r° |0 s (x)| 2 ^ +2 rfx = (/i + l)^ 2 ^ ! (l-t 2 )Ut. (3.21) 

J-oo JO 



The proof is concluded. 



□ 

The previous proof shows that the infimum d°(u) is approximated by the action of a soliton escaping 
to infinity. Moreover notice that the minimizing sequence weakly converges to the vanishing function. 

4. Variational Analysis: the 5 vertex 

In this section we discuss the variational properties of the action functional in the general case with 
a < 0. In fact we prove that there exists a* < such that for —N^/uJ < a < a* the action functional 
constrained to the Nehari manifold admits an absolute minimum. The proof of this statement 



is broken into several lemmas. Firstly, in Lemma 4.1 we prove an equivalent formulation of the 



variational problem we are studying. Then, in Lemma 4.2 and Proposition 4.1 we prove that the 
infimum of the constrained action is strictly positive and smaller than the infimum of the Kirchhoff 
action, therefore for a negative enough the infimum is not reached by functions escaping at infinity 



like (3.12), otherwise the two infima would coincide. This is a key ingredient in the proof of the main 



Theorem [TJ where we prove that a minimizing sequence admits subsequences with non trivial weak 
limit. Finally we prove that this limit is the absolute minimum. 
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We recall from the introduction the action functional, given by 

rv 



and the Nehari functional 



l2/i+2 

The Nehari manifold is given by £ £, $ ^ 0, s.t. ^[^j = 0}. It is understood that the above 
functionals are defined on the form domain £. The action restricted to the Nehari manifold will be 
named again reduced action and is defined by 

m = 2^||*||jS£a = S u [9] - (4.1) 

We introduce also the function 

d(u) = ini{S u [*] s.t. $G^,f/0, /„[*] = 0}. (4.2) 

In the following let a < 0. 

We want to prove that for a smaller than a threshold value a* the action S u constrained to the 
Nehari manifold admits an absolute minimum. 

Firstly we give an equivalent formulation of this variational problem. 
Lemma 4.1. The following equality holds 

inf s.t ¥ G £, * ^ 0, I w [tt] = 0} = mf{S[V] s.t. * e 5, * ^ 0, I w [tt] < 0}. (4.3) 
Moreover $ G £ satisfies S[$\ = d(u) and < ^ = d(w) and = 0. 

Proof 

The idea is the following: if a function $ is not on the Nehari manifold and < 0, then by 

multiplication by a suitable scalar, it can pulled on the manifold lowering the reduced action at the 
same time. First notice that by ( 4.1[ ) we have immediately 

inf^ftf] s.t. G £, V ^ 0, J w [tf] = 0} > inf{£[tf] s.t. * G 5, * ^ 0, < 0}, (4.4) 

since and 5 coincide on the Nehari manifold. 
Now take $ G £ such that I w [$] < and define 

\ 11^112^+2 / 

Since I u [®] < then (3 < 1. Moreover by direct computation one has 

0. (4.6) 



Then, using again (4.1), one has 

S u \p&] = S[f3<5>] = (3 W S{<5>] < S[Q) (4.7) 

then 

inf{S w [#] s.t. G £, J w [tf] = 0} ^ inf{£[#] s.t. * G £, J w [tf] < 0} (4.8) 
and identity ( 4.3[ ) has been proved. 
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Notice that if $ minimizes S w on 1^ — then it minimizes also S on I w ^ by ( |4,3[ ) . Suppose now 
that S[$] = d(oj) and < 0. Then defining f3 as above one has that S u [/3$] < S[$\ = d{u) which 
is a contradiction to the definition of d{uj). 

□ 



Lemma 4.2. Assume u > a 2 /N 2 . Then d(u>) is strictly positive. 
Proof 

Firstly we derive an elementary Sobolev inequality for the halfline. Let / G and denote by 

f(k) its Fourier transform. Then, we have 

1 r.t 1 



1/(0)1 < 



\f(k)\dk 



By Cauchy- Schwartz inequality we have 



(4.9) 



1/(0)1' < 7iW\\vm + * 



(4.10) 



If G H 1 (R + ) we can extend it by parity to a function on the line and apply (4.10). In this way we 
finally have 

imi^^'llh^+^Ml^y (4-11) 
Now take $ G 8 then by ( 4.11[ ), we have 

N , N 



i=i i=i v 



Lb 



-||$|| 2 + — 1|$'|| 2 . (4.12) 



Then, using again (4.3) and with a suitable choice of a (it is possible due to the restriction on u) we 
have 

lad 



> 1 



\n 2 + {u-a\a\)\\n 2 -\\ntxi>4n 2 m 



By Sobolev type inequalities we arrive at 



which implies, for a non vanishing function $, 

ll$l 



2^+2 



Since on the Nehari manifold S u and S coincide, we must have d(u) > 0. 



mill (4.i3) 



(4.14) 
(4.15) 

□ 



For any u> > define a* such that —Ny/uj < a* < and 



(1 -t 2 )»dt 



N 

N-Ju 



:i-t 2 )»dt. 



(4.16) 



16 



RICCARDO ADAMI, CLAUDIO CACCIAPUOTI, DOMENICO FINCO, AND DIEGO NOJA 



Notice that a* is uniquely denned since the r.h.s. of (4.16) is a decreasing function of \a*\ whose 
range includes the value J (1 — t^^dt. 



(4.17) 



Proposition 4.1. Let —N^/uj < a < a* . Then 

d{u) < d°{u). 

Proof 



In order to prove (4.17), it is sufficient to exhibit a trial function & £ such that = and 

S u [^] < d°(oj). Let a be defined as 

1 



arctanh 



/ii/W 



\a\ 



NJu 



Now we consider the symmetric trial function \l/ given by: 

= (j) s (x + a) i = l...N 



(4.18) 



where 4> s is defined by (B.l). By construction ^ G £. Moreover it is straightforward to check that 
$ G V{H) and 

- |* = (4.19) 



Multiplying both sides of (4.19 ) and integrating by parts, one checks that iw^] = that is \I> satisfies 
the constraint. Therefore it is sufficient to evaluate the reduced action using (B.5). One has 



'l-t^dt 



(4.20) 



and the condition S w [^/] < d°(co) amounts to 

TV * 
2 



H_ 



l-t 2 )»dt< / (l-t 2 )»dt 



which holds true by the hypothesis —N^/uJ < a < a* since 



N 
~2 



(i-e)ndt<- 

H z 



;i -t 2 )Mdt 



□ 



Now we can finally prove Theorem [TJ as stated in the introduction. 
Proof 

Let {ty n } be a minimizing sequence, we prove that there exists a subsequence weakly convergent in 
H l . First notice that H 1 !^ 1 1 2^+2 is obviously bounded (see Lemma 4.1). Recall that for $ G £ 

,2 



> — 1|$|| 2 . 



(4.21) 
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Using Q4.21D and I^n) ^ we have 

This implies ||^/ n || ^ c. Using again J^^n] ^ we have also 

WKW 2 ^ ll^nllJS-^ll^nf-al^iWI^c + cll^llll^ll ^ c + c(^\\* n f + e\\V n \A (4.23) 

for any s > 0. Taking e sufficiently small we see that \\^' n \\ is bounded and therefore also H^JJIh 1 is 
bounded. By Banach-Alaoglu theorem there exists a weakly convergent subsequence, which will be 
still denoted by {^/ n }- Let be the weak limit. 

Now we prove that ^ 0. To this aim we preliminarily show that G E, ^ n (0) — > ^oo(O) 
and that J^f^n] 0. Let A 3 ' : Q — > M be a function on the graph defined in the following way: 
Xj(y) = e~ y and \i(y) = for i ^ j. Then by weak convergence and integration by parts we have 

^>(0) = (A j , VJhi ->• (A\ V^hi = * itO o(0). (4.24) 

Since i[)j >n (0) does not depend on j, the first two claims are proved. We prove the last claim by 
contradiction. Assume that 

I u [V n ] ->0 (4.25) 
is false, then there exists a subsequence, still denoted by {^/ n }, such that 

lim J w [* n ] = 7 <0. (4.26) 



Let (3 n be defined according to (4.5) then 



7/i 



2/i 



lim bin [l+ / u[ l n L 1 I I , L f— ) <:i ( 1.27) 

therefore 

lim S[(3 n V n } = lim ffi +2 S[S> n ] < d{u) (4.28) 

n— ¥oo n—>oo 

and /^[/Jn^n] — but this contradicts the assumption that \I/ n is a minimizing sequence. Hence 
0. 

We proceed again by contradiction to prove that 7^ 0. Assume that = and define 

[\\K\\ 2 + ^n\\ 2 ]^ 

hi/ " 



2(i 



Using (4.24), (4.26) and the contradiction hypothesis, one has 

lim p n = lim f 1 + ^"iTltr"' ) - 1 < ! - 29 ) 

\ H v ll2/i+2 / 

Therefore 

lim = lim pl^S^n] = d{u). (4.30) 

n— >oo n— >oo 

On the other hand, by direct computation one has 

J°[p n vp n ] =0 . (4.31) 
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Therefore, by Proposition AA_ and Theorem [3] 

d(u) < Sooiu) < S[p n ^ n ). (4.32) 

Passing to the limit, one obtains 

d(u) < Sooitu) < d(u) (4.33) 
therefore the hypothesis = can not hold. 

Now we shall prove that i^f^oo] ^ 0. We recall, see [12], Brezis and Lieb's lemma: if f n converges 
weakly to foo in L p , 1 < p < oo, then 



In our case, this implies that 



li/»lls-ii/»-/ooii;-n/ooii5->o. 



and, applying (4.34) both to and that 

~ - ^oo] " Woo] -)• 0. 

Suppose that i^oc] > 0. Then, by fl4~25) and {OfJ, 

lim J w [tf„ - = lim J w [*„] - J^oc] = < 0. 

Choose n such that i^^n — < for n > n. Then by definition of d(u) we have 

d(u) < S[^ n - ^oo], n > n. 



On the other hand, since ^ 0, by (4.35) one has 



lim S[V n - ^ = lim S[* n ] - Sptfoo) = d(u) - S^^] < d(u) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



and this contradicts (4.38); so it must be i^^oo] ^ 0. 

By definition d(u) ^ S'^oo]. On the other hand, by the lower semicontinuity of the norm under 
weak convergence we have 



1 J 2(/i+l)" °° ll2 ^ +2 ^^cx 



s II *n II %ll = d(u) 

2(/i + l)" 11211+2 v ; 



which implies 



SI*,*] = d(u) 



(4.40) 



(4.41) 



and so is an absolute minimum of S u constrained to Nehari manifold. 



□ 
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5. Stationary States 

In this section we explicitly compute the stationary states of S u and of S® and identify the 
minimum of the action. Some of the results in this and subsequent section were announced in |4]. 
We denote by [s] the integer part of s. 

Theorem 4 (Stationary states of S u ). 

Let a < and cu > f^; then S w has [(N — l)/2] + 1 critical points ^ u> j, with j = 0, • • • , [— jp' 1 ], given, 
up to permutations of the edges, by: 

[0 a (a; + aj) z=j + l,...,iV 

a, = — = arctanh ( — —) . (5.2) 

Moreover, for —N^/uJ < a < a* the function \& w o is the ground state. 



Proof 

A regularity argument shows that a constrained critical point of the action S u is in fact an element 
of the domain of the operator H . We sketch the standard proof. Any such non vanishing critical 
point * satisfies S'(&) = 0, i.e. 

S^(*)$ = 0, V$ G 5. (5.3) 



Applying (5.3) first to <3>, then to H = — z$, and summing the two expressions, we find 

J B(vI/,$)--^-(|vI/| 2 ^ ) $) +w (vI/ ) $) = 0, (5.4) 



where B is the bilinear form associated to the quadratic form E hn . So, from (5.4) the following 
estimate holds 

< CW||$||, V$ G £. (5.5) 

Notice that, choosing $ among the functions vanishing in a neighborhood of zero, we conclude from 
(5.5) Riesz theorem and definition of weak derivative that every ipi G H 2 (R + ). Thus, for a generic 
<J> G £ an integration by parts gives 



N / N 



2 =-^(^ ,/ ,0 4 ) L2(K+) -0i(O) o;Vi(0)-^V/(0) . (5.6) 



8=1 \ 8=1 



So, from (5.5) and ( |5.6 ), we conclude that \I/ belongs to the domain T)(H). Moreover, the function 

Hm - 2^\^\ 2 ^ + belongs to L 2 (R + ). 

Therefore S'^ = is equivalent to the following equation 

ff* w -|* w | 2 "* w = -w* u w>0. (5.7) 

Notice that H acts locally as the Laplacian, thus on every edge we must seek L 2 (M + )-solutions to 
the equation 

-0" - |0| 2m = -u<j> co>0. 
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The most general L 2 (M 
a G C, |er| = 1 and y G 



)-solution is <p(x) = <r(f) s (x — y) 
:. Therefore the components {^> L 



i i 

a [(// + 2 ^ sech^ (fj,y/uj(x — y)) where 

j of a critical point ^ are given by 



{x) = o % (\){x - y{) 



In order to have a solution of (5.7) it is sufficient to impose boundary conditions (1.2) such that 
ty u G T>(H). The continuity condition in (1.2) implies 0\ = . . . = and y^ = e^a with e$ = ±1 and 
a > 0. We can omit the dependence on o without losing generality. Referring to the bell shape of 
the function <p s , we say that in the i-th edge: there is a bump if y^ > 0, that is, if £j = +1; there is 
a tail if yi < 0, that is, if £j = — 1. Now we determine £j and The second boundary condition in 



(1.2) rewrites as 



tanh(/x 



N 



Co' 



(5i 



Equation (5.8) gives as a first constraint that Xlili £ « must have the same sign of a. That is the 
critical point must have more tails than bumps. For every such a configuration, or equivalently a 



choice of the set {£i}, condition (5.8) fixes uniquely a. We choose to index the solutions by the 
number j of bumps. Correspondingly one obtains a unique solution to ( |5.8[ ) which we call aj. In 
this way we arrive at (5.1) and ( |5.2 ). For instance, if N = 3 then there are two stationary states, a 
three-tail state and a two-tail/one-bump state. They are shown in figure [TJ 





Figure 1. Stationary states for iV = 3, a < 



To summarize, solutions to (5.7) are given by \P 



:u))3 with j = 0, 

■ ■ 2 

Notice that (5.8) admits solutions iff the lower bound ^ < u 



...,[(#- l)/2]. 

holds true. We can explain this fact 
for a < noticing that (5.7), for small ^ and neglecting nonlinearity, is the eigenvalue equation for 
the linear part of the Hamiltonian corresponding to energy E = —u; taking into account the known 

2 

fact that the linear graph Hamiltonian H has the ground state energy —^2, the lower bound means 
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that the nonlinear standing waves bifurcate from the vanishing wavefunction at the ground state 
energy of the linear problem. 

Now we prove that \l/ Wi o is the ground state. Notice also that ^,0 is uniquely defined since it is 
invariant under permutations of the edges. We know that for —N^JuJ < a < a* a minimum of 
exists and therefore it is a critical point. It is sufficient to prove that S^f^^o] < S^l^^j] for j ^ 0. 
In fact we prove a stronger statement, that is, if ^ j ^ [(N — l)/2] — 1 then 



(5.9) 



Using (B.5), equation (5.9) is equivalent to 

-1 



J 



(1 -t 2 )»dt + {N -j) 



(N-2j)^n 



[1 - t 2 )ht < 



V + 1) 



Let us define the constant C 



(1 — t 2 )^dt + (iV — j - 1) 



(N-2j-2)^J 



[l-t 2 )»dt. (5.10) 



(JV_2j-2)v^J 



c 



'l-t 2 )-»dt. 



It is convenient to rewrite the l.h.s. of (5.10) as 

-1 



J 



(1 -t 2 )»dt + (N-j) 



(N-2j)^n 



{l-t 2 )»dt = 2jC + (N -2j) 



(JV-2 3 ')\A5 



;i -t 2 y»dt 



(N-2j)V^ 



|q| 

NC-(N-2j) / (l-t 2 )»dt. (5.11) 

Jo 



Repeating the same manipulations for the r.h.s., we see that (5.11) is equivalent to 

I a I I a I 

(N - 2j - 2) / (1 - t 2 )»dt < (N — 2j) / (1 - t 2 )»dt. 

Jo Jo 

With a straightforward change of variables the last inequality becomes 



(5.12) 



t 



N-2J-2 



dt < 



1 - 



t 



N-2j 



dt 



(5.13) 



which is manifestly true. 



□ 



Remark 5.1. For N > 2 and j > there exist excited states, but only for parameters 00 > ^ N 1 2 j) 2 • 

2 

So the picture is that for fixed a and increasing lo firstly the branch of ground state is born at lo > 
and then for sufficiently high u the branches of higher excited states appear. 



22 



RICCARDO ADAMI, CLAUDIO CACCIAPUOTI, DOMENICO FINCO, AND DIEGO NOJA 



Remark 5.2. Even if in the present section we considered the case a ^ 0, notice that the analysis 
of the previous theorem can be repeated also for a > and one would find that the critical points 
are given again by (5.1) and (5.2) with j = [N/2 + 1],...,N. This means that for a repulsive 5 



interaction at the vertex the stationary states have more bumps than tails. 
We end this section with the characterization of stationary points of 5° 

Theorem 5 (Critical Points of S®). Let u > 0. If N is odd, then there is a unique critical point 
of S® given by 

= & i = l,...,N (5.14) 

If N is even then S® has a one parameter family of critical points given by: 




l,...,N/2 
N/2 + 1, . . 



N 



a G 



(5.15) 



Proof 



Repeating the argument in the proof of Theorem |4j we have to find the solutions of 

N 

t&nh(ny/uja) = 0. 

i=l 

If TV is odd, then there is a unique solution given by a = which corresponds to the critical point 



(5.14), and such a solution can be described as composed by N half solitons glued at the vertex. 

,...,N/2, 



On the contrary, if N is even then there are infinite solutions: a G 

Si ■■ 



Si = +1 for i 



■1 for N/2 + 1, . . . , N gives a solution to (5.15) which corresponds to \I/°' a . 



□ 



Remark 5.3. If N is even, then the graph can be considered as a set of N/2 copies of the real line. 
With a Kirchhoff boundary condition, one has continuity and derivability of the wavefunction at the 
vertex, and the above solutions ^°j a can be interpreted as N/2 identical solitary waves on each real 
line translated by a quantity a. 

Remark 5.4. In the case N = 3 and for a cubic nonlinearity it has been proved in j3] that the energy 
E at constant mass M is not minimized on , which turns out to be a saddle point. In fact, the 
constrained energy is bounded from below but it has not an absolute minimum. We conjecture that 
the same phenomenon happens here for the action. 

6. Stability of Ground States 

In Section 4 we showed the existence of a profile ^uj,o (denoted there by \IZoo) which minimizes the 
action for the star graph with attractive delta boundary conditions at the vertex. This minimizer 
is the ground state of the problem if the strength a of the point interaction at the vertex is sufficiently 
large. In Section 5 we provided the explicit expression of stationary states f WJ , and in particular 
of the ground state ^l/^o- In correspondence to the ground state (and to every stationary state) one 
has a standing wave of the form ^ e tLUt which solves the NLS on the graph. In this section we 
study the stability of such a standing wave. Being a time-dependent solution and not an equilibrium 
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point of the autonomous equation (1.4), stability has to be intended as orbital stability. This means 
Lyapunov stability up to symmetries of the equation, which in this case are related to the gauge 
£7(1) invariance of the Hamiltonian of the problem. To be precise, we recall that the orbit of *f> w is 
defined as ^(tfj = {e ie ^ u , 9 G K}. 

The state is orbitally stable if for every e > there exists 5 > such that 
d(tf(0),tf(tf w )) < 6 d(tf(t), < e W>0 



where \I/(t) is the solution to (2.1) with initial data \l/o 



d(V, = inf l|*-<%, 

and the norm || • ||g is the energy norm, given in our case by H 1 norm in S. 
A stationary state is unstable if it is not stable. 

Orbital stability of solitary (not necessarily standing) solutions to nonlinear Schrodinger equations 
is a well developed subject, studied in several classical papers. Two main techniques have been 
developed to establish orbital stability of solitary waves: the method of Cazenave and Lions based 
on Concentration Compactness ([TJJ [15]), and the method of constrained linearization pioneered by 
Benjamin in the case of KdV equation and studied more sistematically by Weinstein and Grillakis- 
Shatah-Strauss [351 |36j [191 120] • In a forthcoming paper we will treat minimization of energy at 
constant mass through a suitable adaptation of concentration compactness method to the case of 
star graphs ([I]), while here we refer to the Weinstein and Grillakis-Shatah-Strauss method which 
is especially suited for treating stability of equilibria of Hamiltonian systems with symmetry. Some 
preparation is needed to cast our problem in this framework. As in the scalar case, the NLS on a 
graph turns out to be a Hamiltonian system on the real Hilbert space of the couples of real and 
imaginary part of the wavefunction. We pose \I/ = U + i V = (U,V), where U = («!,••• ,un) t 
and V = («!,■•■ ,v N ) T . So we identify L 2 {Q) = L 2 (£,C) with L 2 (Q,R) © L 2 (Q,R) := L|(£) . 
Analogously one can define the spaces L^{Q). Correspondingly, L 2 (Q) can be given the structure of 
a real Hilbert space taking as its scalar product the real part of the usual complex one: 

((U u V.f, (U 2 , V 2 ) T ) Ll{g) = Re^, m 2 ) L2{g) . (6.1) 

Furthermore, L 2 (Q) is also a symplectic manifold when endowed with the symplectic form (coinciding 
with the imaginary part of the complex scalar product) 

- r 

n((U u V{), (U 2 , V 2 )) = Im^x, ^ 2 ) L 2 (g) = Y] / (0>2)i(«i)< - (vi)i(u 2 )i)dx. 

The same symplectic structure is inherited by the energy space S. Moreover, multiplication by the 
imaginary unit i is equivalent to acting by the matrix — J G Mat(R, 2N x 2N), where 



J 



/ 
-I 



and the blocks and / are the zero and unit matrices in Mat(R, iV x N). 

Note that if \I/ G V(H), then the real vectors U and V satisfy the same boundary conditions as 
we will say, with a slight abuse, that they belong to T>(H). With these premises, the nonlinear 
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Schrodinger equation for \l/ is equivalent to the canonical system 

d_ fu" 

dt [y, 

where the Hamiltonian E becomes 



JE'[U,V\. 



mV) = -||(?7',V)||5 fi(W - 



^ll(^^)||^ 2(g) + f(K(o)P + K 



(6.2) 



(6.3) 



or explicitly 



N 



k=l 



(Kl + kl ) dx 



a .. 



AT r 



2/x + 2 



E 

k=l 



(\u k \ 2 + \v k \ 2 r +1 dx 



2 + l^i(0)| 2 ) 
= E[u,v) , 



(6.4) 



and the derivative E' is given by 



E'((U,V))[(H,Z)) = -{E((U,V)+e(H,Z))} €=0 . 

Linearization of the Hamiltonian system ( |6.2 ) around the stationary state is achieved by substituting 

(*t)j = ({* u ,o)j + h j + iz j )e iut (6.5) 
and neglecting higher order terms than linear in (6.2). The real vector functions H and Z satisfy 

d fH s 



dt 



JC 



where £ is the unique s.a. operator associated to the symmetric and lower bounded quadratic form 
£"($ Ui o), i.e. the second derivative of the energy at the ground state. Indeed, the second derivative 
is defined through the formula 

C((H 1 ,Z 1 ),(H 2 ,Z 2 )) = E"(^ u}fi )((H 1 ,Z 1 ),(H 2 ,Z 2 )) = J^{E(* u + e(H 1 ,Z 1 )+\(H 2 ,Z 2 )} f ^,x = * • 

An easy computation shows that C = diag(£_,£ + ) and the matrix operators £_ and C + are given 
by (here the summation convention is used) 

d 2 



dx 2 
dx 2 



+ UJ 



,0)k\ 



2M kk 



+ cj-(2 f i+l)\(* ufi ) k \ 2 ' i )8 i , k . 



The operators C- and C + act on the real vector functions H and Z belonging to T>(H) = T>(C±). 
Notice that to simplify notation from here on we suppress the dependence of operators C± on the 
ground state ^l/^o- 

Precise conditions to have orbital stability (and instability) for general Hamiltonian systems and in 
particular for systems of NLS equations, are given in the already quoted papers of Weinstein and 
Grillakis-Shatah-Strauss. They can be reduced to the validity of three conditions, called Assumptions 
I, II and III in [19J and j20] and the verification of a further convexity condition on the function 
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d(oS) = S w (^ w ) introduced in Section 4 and called in the physical literature the Vakitov-Kolokolov 
condition. 

Assumption I is the well posedness, proved in Section 2. Assumption II is the existence of a regular 
branch of standing solutions of the stationary equation, proved for our model in Section 4 and 5, where 
the regular family of standing waves (J^,+oo) 3 u h-> ^ Wj0 is explicitly constructed. Assumption 
III concerns spectral properties of linearization £"'(^ W)0 ) = (£_,£+) around the ground state. The 
spectral conditions are stated and proved in the following proposition. 

Proposition 6.1. The operators and C + are self adjoint. Moreover: 
i\) kerC + = {^,0} an d the rest of the spectrum is positive; 
i 2 ) ker £_ = {0}'; 

is) n(C ) = 1, where n(A) is the number of negative eigenvalues of the operator A, i.e. its Morse 
index. 



Proof 

We begin to remark that operators £_ and £+ are selfadjoint on V(C±) = V(H), due to the fact 
that the components of the ground state (^/^,o)fc are continuous and strongly decaying at infinity, 
and as such they constitute in the matrix operators C± a relatively compact perturbation of H + u. 
For the same reason, by Weyl's theorem the absolutely continuous spectrum of £_ and C + coincides 
with the essential spectrum of H + uj, i.e. [u, +00), and the discrete spectrum is composed at most of 
a finite number of eigenvalues. Let us consider the kernel of C + . This surely contains ^ Wt o- Indeed, 
the equation £+^^,0 = coincides with the stationary equation satisfied by ^ Wi q. Let us show that 
there are not other elements in the kernel. An integration by parts allows to rewrite the quadratic 
form of C + , for any element V G T>(C+), as follows 

N r +oc d 7) 

N , ( 6 - 6 ) 

+ 2^ wKi ) - 



k=l 



(**,o)fc(0) 



and the last term is vanishing due to the 5 boundary conditions at the vertex, continuity and 

Ef =1 ^(o) = OTl (o). 

So (C+v,v) L 2(g) > for every v G T>(£ + ) not coinciding with Nl/^o, which is the only eigenvector 
with eigenvalue of the operator C + . This proves statement i\. 

Concerning statement 12, it is sufficient to consider the equation C-U = 0. This is written, in 
components, as 

d 2 

~ d^ Uk + UJUk ~ ( 2 A t + 1 )l(^.o)fc| 2M Mfc = k = l---N, (6.7) 

where «i(0) = 7/2(0) = ■■■ = un(0) and 52 fc=1 u' k {Q) = au%(0) due to boundary conditions. The 
general theory of second order differential equations gives for the previous equation a solution which 
is a linear combination of asymptotically exponential fundamental solutions; for x — > +00 only one 
of them is in L 2 . Now notice that the function Uk ■ (0, +00) 1— >■ K s.t. Uk(x) = ■^(fy t j,o)k(x), k = 



1, • • • ,N satisfies equation (6.7) in (0, +00) and it is decaying at infinity. So every solution in 
u G L 2 (0, +00) of (6.7) is a multiple of such a function: Uk(x) = CkUk(x),Ck G 5L To conclude, a 
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direct calculation shows that the only real constants compatible with the boundary conditions on u 
is Cfc = 0, k — 1, • • • , N. This proves statement i<i . 

Let us consider finally the Morse index of the £_ operator. It is immediate that n(£ ) ^ 1 . Indeed, 
let us consider the quadratic form for £_ evaluated on the ground state: 

(£-* w ,o,* w ,o)j2(s) = (£ + ^,o,^,o)Ll(g)-2/i(|^,o| 2At ^,o,^,o)^(g) = 0-2||* w>0 ||^ < . (6.8) 

So the s.a. linear operator £_ has a negative vector, so it surely admits at least negative eigenvalue. 
Let us prove that it has a single negative eigenvalue only. This is a consequence of the variational 
properties of ^,0- in fact is a minimum point of the action S w on the codimension one constraint 
I w = 0. This minimization property entails that S"(^ ult0 ) is positive definite on the tangent space 
at fu )0 of the constraint manifold. Being the constraint a manifold of codimension one, S'^^^fi) 
admits at most one negative eigenvalue and the same is true for its only possibly negative diagonal 
component £_. See Appendix B in [T5] for the detailed argument. 



□ 



The last property needed to show orbital stability of the ground state is the so called slope 
condition, or Vakhitov-Kolokolov condition. This coincides with the convexity of the function d(u), 
or more explicitly it means that on the branch of stationary solutions {\I/ W) o} parametrized by u, one 
has 

d 2 „ _ , d 2 , d 



d"\ 



iE(* u )+u>M(* u )) 



I* 



w.O 



> 



du A duj A du 

In fact, a direct calculation making use of the formulas in the appendix (and which is possible in this 
model due to the explicitly known form of Nl/^o), gives 



d 

duo 



Ldfl 



with C = C(N,n,u) 



N 



C 



,1 



—l 



dt + 



a 



2N 



1 - 



a 



N 2 oj J 



(6.9) 



Now, the r.h.s of (6.9) is positive thanks to the lower bound onw > ^g-. The Vakhitov-Kolokolov con- 
dition with the spectral properties proved in proposition |6.1 thanks to the Weinstein and Grillakis- 
Shatah-Strauss theory constitute the proof of Theorem [2] stated in the introduction. 

Remark 6.1. Note the following facts. 

The theorem gives orbital stability of the ground state also for the critical nonlinearity \x = 2 . 
From formula (6.9) it follows that for supercritical nonlinearities fi > 2 there exists u* such that 

2 

is orbitally stable for to E (f^,ui*) . In [19J it is shown that if Assumptions I, II, III are satisfied and 
d"(u>) < 0, then the standing wave corresponding to u> is orbitally unstable. Again from formula (6.9) 
we see that for u > u* the ground state ty u>0 is orbitally unstable. The case u = u* where d"(u>) = 
is undecided. 



Appendix A. Rearrangements 

For a given function $ : Q — > C we introduce the rearranged function <3>* : Q — > R. The function $* 
is positive, symmetric, non increasing and is constructed in a way such that it is equimisurable w.r.t. 
$, that is, the level sets of |$| and $* have the same measure. This is sufficient to prove that all the 
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L P (Q) norms are conserved by the rearrangement. The comparison of the kinetic energy of $ and 
$* is more delicate. On the real line the Polya-Szego inequality shows that the kinetic energy does 
not increase. This is no longer true for a star graph where a constant N/2 appears, see Theorem [6] 
below. 

Given $ : Q — > C, we introduce A(s) and /i(s) defined by 

K') = IOT > 41 M*) = IOT>«}I- 

Now we define the symmetric rearrangement of $. 
Definition A.l. Define g : R + — > R + as 

g(t) = sup{s| A(s) > Nt} , 

i/ien we put $* = (0*, 4>* N ) with 

<j)*{x)=g{x) j = l,---,N. 

The main properties of $* are the following: 

Proposition A.l. Let $ 6 L P (Q) . The symmetric rearrangement $* zs positive, symmetric and 
non increasing. Moreover = 

Proof 

By construction $* is positive and symmetric. Since A is non increasing, $* is non increasing too. 
Now we prove the invariance of the LP norm. First we prove that if ^ Nt ^ /i(0) then 

V(g(t)) ^Nt^ X(g(t)). (A.l) 

Let s' < g(t) then, by definition of g, A(s') > Nt. Since the latter estimate holds for every such s', 



taking the supremum we have X(g(t)) ^ Nt which is the second half of (A.l). 

Now choose s' > Nt. Then /z(s') ^ A(s') ^ Nt and taking the infimum over s' the proof of (A.l) is 
complete. 

A key property of the symmetric rearrangement is the equimisurability, that is: 

\{$*>s}\ = \{m>s}\. (A.2) 
Notice that the set {|$| ^ t} can be rephrased as the union of disjoint sets 

{|$| >t} = {\(h\ >t}\J...\J{\<j> N \ >t}, 
where {\<fij\ ^ t} is to be understood as a subset of the j-th edge. Therefore we have 

TV 

\{m>t}\ = Y,\m>t}\. 

i=l 

Fix s and define 

t = swp{t\g(t) ^ s} . 



If g{t) > s then Nt ^ X(g(t)) ^ A(s) by (A.l). Taking the supremum over t we get Nto ^ A(s) 



Assume by absurd that Nt < A(s). Take t such that Nto < Nt < A(s) then g(t) ^ s and the 



contradiction with the definition of to is reached. Since Nt = |{$* ^ s}\ equality (A.2) is proved. 
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By the layer cake representation (see [27] ) we can prove the invariance of L p norm under rearrange- 
ments. We start from 

\m\l = J2 \My)\ p dy = P s^- l \m^s}\ds 

• =1 JM.+ JO 



and using (A. 2) we obtain 



+oo 



H P P = P > s}\ ds = ||$*||£ 



which is the desired identity. 



□ 



Notice that if {|$| = t} has non-zero measure for some t then A has a jump, while if |$| has a jump 
then A has a constant part. Moreover if A has a constant part then g has a jump and if A has a jump 
then g has a constant part. 

Notice also that if |$| is continuous and {|$| = t} has zero measure for t then g is the inverse of A 
up to scaling. 

Now we turn our attention to the Polya-Szego inequality and prove it by elementary methods. We 
mainly follow [2T] while some technical results are taken from |22j . 

From now on, we restrict ourselves to real and positive $. Later we shall extend the result to the 
general case. First we gather some preliminary results in Lemmas A.l and A. 2 Then we establish 



the required estimate for a class of regular functions (Lemma A. 3 and Lemma A.2). Finally we give 



the main theorem the proof of which relies on a careful decomposition of the kinetic energy and on 
a limiting argument. 

Moreover, since functions in £ are continuous, for our purposes in the following we always assume 
that $ is continuous without losing generality. 



Lemma A.l. 

$* e u\g). 



Assume that $ : Q — > R + is continuous and <£> G L P (Q) then is continuous and 



Proof 



Due to proposition A.l 



we have to discuss the continuity part only. Since $ is continuous, A is 
strictly decreasing and may have at most a countable number of discontinuity of first kind. There- 
fore A is locally continuous away from discontinuities and g is locally (up to an irrelevant scaling) the 
inverse function. Then g is locally continuous by the inverse function theorem. At the points where 

that g has a constant part joining the 



A has a discontinuity it is easy to check, using definition A 
non constant branches. So that g is globally continuous. See 



for more details. 



□ 



Reasoning as in [2T] we get the following: 
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Lemma A. 2. Let $ n , $ : Q -)■ IR + and $ n , $ G L P (Q). Suppose that ||$ n — <&|| p ->■ 0, i/ien 

^ liminf A n (s) ^ limsup A„(s) ^ A(s) , 
» n 

0(0) < lirninf<? n (0). 

n 

We introduce the following class of regular functions. 

Definition A. 2. Let VL be the set of functions $ : Q — > 1R + such that: $ is continuous, compactly 
supported and, for any j = 1, ...,N, there exists a finite number of compact intervals J Jjn such that 
supp = M J^n and $ restricted to is affine. 

This class of piecewise linear functions is dense in S. 

Lemma A. 3. Let $ G VL. There exist two open sets 0\ and O2 such that $* is constant on 0\ and 
$* is differentiable on O2 with |$*'| > 0. Moreover, Q \ (0\ U O2) consists of finitely many points. 

Proof 

Let = a < ai < ... < a m be the values assumed by $ at the boundary of all Ij >n . If the set 
Ui{$ = ai} has positive measure then $' = a.e. on such a set. In the same way $*' = a.e. on 
the set Ui{$* = a;}. We define Oi to be the interior part of U;{$* = a*} and 2 = Q \ Ui{$* = a*}. 
By construction £ \ (Oi U O2) consists of finitely many points. 

We have to show that $* is differentiable on O2 and |$*'| > 0. We introduce the notation Di = 
{ai-i < $ < ai} and D| = {aj_i < $* < aj}. Each set .D; is decomposed first into the components 
on each edge, that is, Di = UjDj. Then we further decompose each D\ into a finite union of open 
intervals D\ k = 0, 1, . . . , n — n(i,j), such that $ restricted to D\ k is affine and non constant (see 



2|). Let us fix s such that aj_i < s < a^. Then the equation $ = s has a solution y J ik (s) G -D 5 



Figure 

for each k = 1, . . . , n. We enumerate the sets D 3 ik 'm.k for i,j fixed in an increasing way w.r.t to the 

distance from the vertex such that y J il < y\ 2 < . . . < y\ n . We put n(i,j) = if D\ = and no y{ k 
is defined for that values of i and j. 

We introduce also D* defined as the projection of D* on the first edge. Let y*(s) G D* be the solution 
to g — s. We can express the measure of the level set {$ > s} by means of the y\ k . One has 

JV N n 

|{$ >S}\=J2 Vin - VU-I + vi,n-2 ~ 24-3 • =£ D" 1 )"^ 
j=l j=l k=l 



Therefore by (A. 2) one has 

N n 

|{$* > s}\ = N\{g> s}\ = Ny* = Y, £(-l) n+ %- ( A ' 3 ) 

3=1 k=i 

On each set D\ k the derivative $' does not vanish by construction which implies that the function 
y J ik (s) are differentiable w.r.t s and 



dy: 



-1 



$' = I I on D 

as 



j 



,k 
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Figure 2. Definition and enumeration of D 3 ik 



by the inverse function theorem. The function y*(s) is differentiable by equation (A.3). It is also 



invertible since s is away from values of $ corresponding to level sets with non zero measure. Therefore 
for such values of s the function g is strictly decreasing. By the inverse function theorem y* is 
invertible and 

^0. 



9 



<kf_ 
ds 



Let L be the Lipschitz constant of Adapting the reasoning in 
estimate: 



□ 

one can prove the following 



dy* 



ds 



This estimate provides an upper bound on g' . Notice that the proof actually shows that $* G VL 
since it says that y*' and therefore g' is locally constant on each D*. If $ is smooth, say C h , then 
the same property holds on O2 for $* by the inverse function theorem. 



Proposition A. 2. Let $ G VL. Then, 



N 

|$*'|| < _ 11$' I 
1 2 



(A.4) 



Proof 

We shall use the notation of the previous lemma. First we consider the r.h.s. of (A.4). We can 
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restrict the integral to the region where $ is not constant and change the integration variable. 



1$ 



/||2 



m N n 

EEE 

i=i j=i k=i 



m N n 



1$ 



/|2 



^ i.k 



EEE 

1=1 j = l k=l Ja i-1 



dy{ 



-2 



ds 



dy\ 



ds 



ds. 



We can repeat the same operation for the l.h.s. of (|A.4|) 

\\&*'\\ 2 = N 



\9 



i\2 



[ ai (dy*'\ 


-2 


dy* 


V d8 J 




ds 



ds. 



Now the conclusion follows by using (A. 3) and the convexity properties of the square function. 
First notice that 



dy* 


1 


ds 


N 



N n 
3=1 k=l 



+k d yU 



ds 



I N n 

EE 



N 



3=1 k=l 



dyj 



ds 



dy 3 

The restriction of $ to an edge has a seesaw behavior and —r^- has an alternating behavior in k. 



Therefore in order to prove (A. 4) it is sufficient to show that 

\ 2 



N 2 



N n 

EE 

3=1 k=l 



( 



d vl,k 



ds 



ds 



( 



N n 

EE 

j=l k=l 



N 



ds 



N n 

EE 

1=1 h=l 



dyL 



ds 



which is equivalent to 



N n 



iEE 

3=1 k=l 



\ 



d vl,k 





ds 




1^1=1 l^h=l 


d v\, h 

ds 



EE — 

dy ] 



3=1 k=l 



ds 



(A.5) 



By the convexity of the square function, inequality (A.5) holds true if 



\j=i fc=i, 



> 1. 



(A.6) 



Notice that X]j=i J2k=i represents the number of solutions to the equation $ = s for aj_i < s < a, 
on the whole graph. Since $ is continuous and compactly supported, there are always at least two 



solutions and then (A.6) holds true. 



□ 



Theorem 6 (Polya - Szego inequality). 

Assume that <3> e S then $* G £ and 



2 11 1 



(A.7) 
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Proof 

Due to Proposition A.l $* is symmetric, and then continuous at the vertex. So that it is sufficient 
to prove (A.7). Let $ G £ be positive and take $ n G VL such that $ n — > $ in H l (Q). Take also 
a positive test function x 6 Co°(<?)- Moreover, let = a < a\ < a 2 < . . . be the values such 
that {g = a>i} has strictly positive measure. Notice that g restricted to D* = {a i _ 1 < g < a*} is 
monotone and invertible by Lemma |A.l Monotonicity of g also implies that its derivative exists 
almost everywhere and is in Ll oc (M> + ). Then the following inequalities hold: 



_ ' poo 

Xj(y)g'(y)dy 

3=1 J ° 

N 

-z2z2 I xAy)g\y)dy 

j=l "* Di 

Yl / Xj(y*(s))ds 



i=i 



N /-ffn(O) 

. liminf V / X j(y* n { s )) ds 

3=1 J ° 

= liminf 

= liminf \(x,K')\ 

n 

^ liminf ||$*'|| 

n 

N 

< — liminf ||x||||$J| 

Z n 



N, 



x\\W\\- 



(A.8) 
(A.9) 
(A.10) 

(A.ll) 

(A.12) 
(A.13) 
(A.14) 

(A.15) 

(A.16) 



The chain of inequalities stands for the following reasons. In (A.8) we explicitly wrote the r.h.s.. 



In (A.9) we have restricted the integral to the regions where g is not constant. In (A. 10) we have 
changed variable of integration, the new one being y*(s) defined as before. This is possible due to 



the restriction made in (A.9). In (A.ll) we have used Fatou's Lemma and Lemma A. 2 In (A.12) 



we have changed back the integration variable and in (A.13) we just changed a sign. In (A.14) we 
have used Cauchy-Schwarz inequality. In (A.15) we have used lemma A. 2 In (A.16) we have used 
the convergence hypothesis. Estimate (A.7) for a positive function / follows from equation (A.16) 
by Riesz Theorem. 



Now we extend the inequality to the general case. First notice that Proposition |A.1| and Proposition 
A.2 both extend to the non-positive case and to the complex valued case. A careful inspection of 
the argument used above, shows that for positive x 



Co°((?) it is still valid until inequality (A.14) 



Then to conclude the proof we have to extend Proposition A.2 to complex valued functions. 
For the real valued case, the extension is trivial. If $ is piecewise linear then |$| G VL and has the 
same rearrangement of $. Notice also that |$'| = ||$|'| almost everywhere. Therefore since (A. 4) 
holds for |<E>|, then it holds also for $. 
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Now for the complex valued case, to define the class of approximating functions we set $ = e l& F with 



6, F e VL, and where the product must be understood componentwise e F = (e 1 f\, e N f 
This set is still dense in S . Again notice that $* = F*. We have also 

|$'|2 _ |Q'_p|2 + |_p/|2_ 

Therefore we can write 

||$*'i|2 = \\F*'\\ 2 < IIFII 2 < 116'FII 2 + IIF'11 2 = ||$'|| 2 



N 



which proves equation (A. 7) in the general case. 

□ 

Remark A.l. The same argument used to prove Theorem^ can be used for the W 1,p norm since 
z ^ is convex for p ^ 1. 

Remark A. 2. The constant N 2 /4 is optimal. For instance take $ such that 

' x < x < 1 
l ) l (y)={2-x l^x^2 (pi = 0fori^l 

x > 2 



TTien g can 6e easz/?/ computed by using ( |A.3[ ) . One /ias: 
From which 



I- fx 0<n<| 
o£ 



7V 2 

ll$ / ll 2 = 2 ||$*'|| 2 = — 
ii ii ii ii 2 

Appendix B. Useful identities 

In this section we recall some useful identities that will be used several times in the paper. We 
label the soliton profile on the real line as 

4> s {x) = [(/i + 1)oj]^ sech^ {p:\fujx). (B.l) 

It satisfies the equation 

- <f>: - \<p s \ 2 y s = -u<p s . (B.2) 
Moreover, multiplying by <p s and integrating one checks that 

II0'JIl2(r) - H0«llK+a(w + w ll^lli 2 (K) = °- ( B - 3 ) 



Starting from definition (B.l) and changing variables in the integrals, one obtains the following 
formulas: 



i , 



x + 0\ 2 dx = (/i + 1) V ^ / _(l-t 2 )^ 1 dt (B.4) 

t 1 Jtanh(^t v / w) 

°°\Mx + 0\ 2>M+2 dx= ^ + 1)1+ " co^ f (l-t 2 )Ut. (B.5) 

Z 1 Jtanh(^ v / ZJ) 
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